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Abstract 

In this paper, the possibility to exploit a high energy beta beam without massive 
detectors is discussed. The radioactive ions are boosted up to very high 7 with the 
neutrino beam pointing towards an instrumented surface located at a moderate 
baseline (e.g. from CERN to the Gran Sasso Laboratories). v e — ► ia, oscillations 
and their CP conjugate are tagged as an excess of horizontal muons produced 
in the rock and tracked by the low-mass instrumented surface installed in one of 
the LNGS experimental halls. We show that the performance of this complex for 
what concerns the determination of the #13 angle of the leptonic mixing matrix is 
comparable with the current low-7 design based on a gigantic water Cherenkov 
at Frejus. 



PACS: 14.60.Pq, 14.60.Lm 



1 Introduction 



Over recent years there has been a marked growth of interest in the development of 
non-conventional neutrino sources. An appealing proposal has been put forward by 
P. Zucchelli |1 in 2002 under the name of "beta beam" . It is based on the production 
of intense beams of /5-unstable heavy ions. The ions are accelerated to a given energy 
and stored in a decay ring with long straight sections pointing towards a far detector. 
Their decays produce a pure and intense v e (z/ e ) beam whose spectrum depends solely 
on the /3-decay kinematics. The advantages of this configuration, the possibility to 
explore subdominant u e — > oscillations at the atmospheric scale with unprecedented 
sensitivity and, hence, to extract the #13 and 5 parameters of the leptonic mixing 
matrix (PMNS [2j) have been discussed by several authors [HI HE]- In particular, it 
has been noted than an European beta beam complex could leverage existing facilities 
at CERN and complement the EURISOL physics program [6 . The latter foresees the 
construction of an intense proton driver for a new generation of radioactive beams. In 
fact, EURISOL is a significant extension of the program presently being carried out 
using the first-generation radioactive ion beam facilities in nuclear physics and nuclear 
astrophysics. It is aimed at increasing the variety of exotic ions produced and their 
yields by orders of magnitude beyond those presently available. Hence, a CERN-based 
beta beam complex would exploit the EURISOL ion source and the CERN PS/SPS 
acceleration complex. Only the dedicated hippodrome-like decay ring should be built 
on purpose. In its most popular configuration v e are produced by 18 Ne ions and v e by 
6 He jl]; the ions are accelerated by the SPS up to 7 ~ 100 ( 18 Ne) and 7 ~ 60 ( 6 He). 
The ratio between the two boost factors is fixed by the equalization of rigidity, i.e. 
the need of accumulating simultaneously into the same ring ions with different Z . The 
corresponding neutrinos emerge with energies below 1 GeV. In order to observe neutrino 
oscillations at the peak of the oscillation probability the detector must be located at 
~130 km from the source, matching, for instance, the distance from CERN to Frejus. 
It has already been noted [5] that low-7 choice is, in principle, quite unfortunate. A 
low-7 beta beam aimed at the observation of — > v e oscillations at the atmospheric 
scale needs a gigantic detector located at L\ ~ 130 km as the proposed 1 Mton water 
Cherenkov at Frejus jjj. The size of the detector must overcome the smallness of the 
cross section at mean v iv) energies of the order of 0.3 (0.2) GeV. A higher energy 
beta beam and a detector located at a farther location L 2 , tuned to operate at the 
peak of the oscillation probability, would provide a flux similar to the low-7 option 
since the neutrino fluxes increase quadratically with the boost factor and decrease as 
L 2 . However, operating at larger 7 show up additional advantages due to the enhanced 
CC cross section, which depends linearly on the neutrino energy. Hence, as a first 
approximation, we expect the sensitivity to the subdominant v e — > channel to grow 
as 7l 2 /7l 1; i-e. as the ratio of the boost factors needed to be at the peak of the 
oscillation probability for the distance L\ and L 2 , respectively. A further increase of 
7 with respect to 7^ 2 would cause a further quadratic rise of the flux compensated 
by a quadratic drop of the oscillation probability ("off-peak" configuration jH]). If the 
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dependence of the oscillation probability on the PMNS parameters were the same in 
the "on-peak" and "off-peak" configuration, this would, anyhow, imply an increase of 
sensitivity due to the further growth of the cross section. Therefore, even under this 
condition, the possibility to run the beta beams in an off-peak configuration could be 
worth being investigated. However, this scenario turns out to be even more attractive 
if we consider the detector technologies that could be exploited to observe v e — > 
oscillations at high z/ M energies. In particular, the v e — > and v e — > z/ M channels 
could be observed as an excess of high energy (> 1-2 GeV) muons from the rock of an 
underground laboratory tagged by an instrumented surface installed into the cavern. 
Since the muon range in the rock grows linearly with the muon energy, the effective mass 
of the rock that contributes to the event rate adds a further linear dependence on 7 so 
that a nearly quadratic increase of the sensitivity due to the higher beta beam energy 
is gained. Clearly, if the background can be kept under control, this configuration 
allows an enormous simplification and reduction of cost with respect to the Mton 
water Cherenkov option, especially if the detector can be installed in pre-existing halls 
as the ones of the Gran Sasso INFN Laboratories. In this paper, we demonstrate 
the feasibility of this design and determine its performance. The structure of the 
subdominant v e — > oscillations for off-peak scenarios and the functional dependence 
of the sensitivity to the parameters of the PMNS on the beta beam energy is derived 
in Sec. 121 The detector concept and the main backgrounds are discussed in Sec. El Its 
sensitivity to the (1-3) sector of the PMNS is computed in Sec. EJ 



2 Oscillations at a high 7 beta beam 

The acceleration of radioactive ions is a prominent technique for nuclear physics studies 
and several facilities have been developed worldwide. The use of short-lived /3-decay 
isotopes transforms these facilities into high intensity sources of pure v e or v e . The 
source has practically no contamination from other flavors and a well defined energy 
spectrum that depends on the kinematic of /3-decay. The annual flux for a far detector 
located at a distance L and aligned with the boost direction of the parent ion is [5]: 



dSdy 



where < y = < 1 - y e , y e = m e /E and 



9(Ve) = ^{y/l- Z/ e 2 (2 - 9y 2 e - 8y 4 e ) + lhy%og 



(2) 



In this formula E represents the electron end-point energy, m e the electron mass, E u 
the energy of the final state neutrino and N@ is the total number of ion decays per year. 
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A beta beam facility based on existing CERN machines has been discussed in [3J E] . 
The protons would be delivered by the Super Proton Linac (SPL) 10J. This driver 
has been studied at CERN in the framework of the neutrino factory [TT] but could be 
an essential part of the EURISOL complex. The SPL would provide 2.2 GeV (kinetic 
energy) protons with an intensity of 2 mA. The targets for ion production would be 
similar to the ones envisioned by EURISOL. In particular, antineutrinos could be 
produced by 6 He decays (a f3~ emitter with E — m e = 3506.7 keV and a 806.7 ms 
half life) from a target consisting either of a water cooled tungsten core or of a liquid 
lead core which works as a proton to neutron converter surrounded by beryllium oxide. 
Neutrinos could result from the + — > 1 + (3 + decay of 18 Ne (E Q — m e = 3423.7 keV and 
half life of 1.672 s); the isotope can be produced by spallation reactions and, in this 
case, protons directly hit a magnesium oxide target. The ions can be further accelerated 
using the EURISOL linac and the PS/SPS complex (see Fig. [TJ and sent to the decay 
ring. In this case the nominal 7 is fixed by the present SPS design. It corresponds to 
7 ~ 60 for 6 He and 7 ~ 100 for 18 Ne. Higher values of 7 can be achieved upgrading 
the SPS with superconducting magnets or making use of the LHC The top rigidity 
available at the LHC could allow for 7 = 2488 ( 6 He) and 7 = 4158 ( 18 Ne) [T2] even if 
in this case the construction of the decay ring would be challenging. In particular, due 
to the high rigidity, it would be unrealistic to build a ring with a curved-over-straight 
section ratio similar to the low-7 option Here, the ratio is 0.94 km/2.5 km (see 
Fig. |U) and the useful fraction of decays is limited by the decays occurring when the 
bunch is located in the return straight section (i.e. the boost has opposite direction 
w.r.t. the far detector). The overall live-time, i.e. the fraction of decays occurring at 
the straight section pointing to Frejus, is 36%. If the size of the straight section is kept 
unchanged (2.5 km) in a high-7 configuration, the live-time is limited by the size of the 
curved parts, which in turns is related to the maximum magnetic field achievable with 
superconducting, radiation-hard magnets. It drops at the level of 20% and 10% for a 
straight section of 2.5 km and a curved section comparable with the size of the SPS 
and the LHC, respectively. The corresponding loss of statistics can be easily and more 
cheaply recovered by instrumenting a larger surface at the far location (e.g. more than 
one experimental hall at LNGS - see below). The actual cost of the decay ring depends 
on the maximum field available and on the possibility to use a significant fraction of 
the accelerating ring also in the storage phase. In the low-7 design 2.9 x 10 18 6 He and 
1.1 x 10 18 18 Ne decays per year are expected. If the LHC were used, some injection 
losses would be expected due to the different optics; these losses could be compensated 
by an increase of the number and of the length of the bunches [5]. 

The number of oscillated v e — > events per year that can be observed at a distance 
L with a detector of M kton is 
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where is the Avogadro's number, 7 is the boost factor of the beta beam complex, 
E is the neutrino energy, e(E) the detector efficiency, a^ c (E) the CC cross section 
at a given energy and P[y e — > Vy) the oscillation probability. The latter depends on 
the baseline L, the energy E and the parameters of the PMNS matrix. In particular, 
in the energy /baseline range of interest, P[y e — > v^) can be expressed as [T5] : 



P(v e -> Vy) ~ sin 2 20 i3 sin 2 23 



sin 2 [(l -A) A] 



+ a sin 2#i3 £ sin 5 sin(A) 
+ a sin 2013 £ cos 8 cos(A) 



1 - A) 2 

sin(iA) sin[(l - A) A] 
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sin(iA) sin[(l - A) A] 
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A 2 

= Oi + 2 (5) + 3 (5) + 4 . (4) 

In this formula A = Am^L/iAE) and the terms contributing to the Jarlskog invariant 
are split into the small parameter sin20i 3 , the 0(1) term £ = cos0i3 sin20i 2 sin20 2 3 
and the CP term sin 6; A = 2y/2G F n e E/Aml 1 with Gp the Fermi coupling constant 
and n e the electron density in matter. Note that the sign of A depends on the sign of 
Am|i which is positive (negative) for normal (inverted) hierarchy of neutrino masses. In 
the following we assume the present best fits for the solar and atmospheric parameters: 
Amli = 7.3 x 10" 5 eV 2 , sin 2 20i 2 = 0.8, \Am 2 31 \ = 2.5 x KT 3 eV 2 , sin 2 2 23 = 1 0. 

Far from the oscillation peak (A <Ci 1, |(1 — A)A\ <C 1) the functional dependence 



5 



of P{y e — > u^) becomes similar to the one of CNGS jH], i.e. 
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and the oscillation probability mainly depends on # 13 and cos 5. Matter effects are 
strongly suppressed and CP asymmetries appear only at the subleading order 2 (6). 
However, even in the highest 7 scenario, the mean neutrino energy remains significantly 
lower than CNGS and the cancellation of the 02(d) term is not complete. Hence, the 
v e flux ( "antineutrino run") 1 still contributes to constrain the (#13, 5) parameter space. 
At fixed baseline L, an increase of 7 implies a quadratic increase of the flux while an 
increase of the mean neutrino energy (E) causes a quadratic decrease of the oscillation 
probability, a linear increase of the cross section (er^ ~ E in the region dominated 
by deep inelastic scattering) and a linear increase of the mean primary muon energy. 
Clearly, 7 and (E) are fully correlated and (E) ~ 7: 



(E) - 

Jo'" 

= 27.^0 — 

z fz^z 2 -y 2 e 

= 2 1 E f(y e ) (6) 

This functional behavior is depicted in Fig. El The upper plot shows the average 
neutrino energy {E) versus 7. The average muon energy for the u^N — > final state 
is also shown. The lower plot represents the number of oscillated events per kton-year 
as a function of 7 assuming 100% conversion probability (Eq.® with P v& _^ u = 1) and 
an isoscalar target. The baseline is fixed at L = 732 km, i.e. at the location of the 
Gran Sasso Laboratories (LNGS). The number of oscillated events (xlO 3 ) per kton-y 
for #13 = 3°, 5 = 0° and normal neutrino hierarchy (Amf 3 > 0) assuming perfect 
detector efficiency is also shown. Due to the assumption e(E) = 1, the linear increase 
of the average muon energy, resulting in a further linear rise of the effective fiducial 
mass, is not exploited. This issue will be discussed in the next section. 



3 The detector concept 

One of the most remarkable features of the beta beams is the clearness of the final state 
to be observed. This advantage is less evident for low-7 configurations (neutrinos in 
the sub-GeV energy range); here, the total cross section is dominated by quasi-elastic 
interactions whose kinematic is obscured by Fermi motion. The latter makes impossible 

lr The antineutrino run is done in parallel with the neutrino one, since it is possible to circulate 6 He 
and 18 Ne ions simultaneously in the decay ring and exploit the time structure of the beam to separate 
the two contributions at the far detector. 
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Figure 2: (Upper plot) Average neutrino and muon energy versus 7 at a baseline of 
732 km. (Lower plot) The continuous line represents the number of oscillated events 
per kton-y assuming 100% conversion probability. The dashed line shows the number 
of events per kton-y for 9 13 = 3°, 5 = 0° and normal neutrino hierarchy (Am 13 > 0). 
Rates are computed for assuming perfect detector efficiency (e(E) = 1). 
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to extract the full spectral information from the final state reconstruction. In this 
energy range, a massive water Cherenkov performing a pure counting experiment could 
be an appropriate detector. At higher neutrino energies, the interactions are mainly 
deep inelastic and the use of denser detectors becomes rather attractive. For a counting 
experiment, a dense tracking detector provides a strong muon/pion and muon/electron 
separation and identifies the neutrino direction through the reconstruction of the muon 
track. Moreover, if the interaction vertex is contained, a kinematic analysis is possible 
and greatly contributes to the determination of the PMNS parameters [SJ. z/ M CC 
interactions occurring into the rock which surrounds the detector can be exploited 
as well. In this case the rock acts as a massive target. Electrons and pions mainly 
interact before reaching the surface of the experimental hall and the tracking capability 
of the detector can be used to veto punch-through hadrons and reconstruct the residual 
muon energy at the entrance of the hall. The main difference with respect to the case of 
fully instrumented volume is that the kinematic analysis is deteriorated by the energy 
loss through the rock but the event rate scales no more as the tracking volume but 
it is proportional to the tracking surface. Hence, an instrumented surface located in 
an existing deep underground hall would imply order-of-magnitude cost reductions 
compared to a massive water Cherenkov. In the following we consider an instrumented 
surface (15x15 m 2 ) installed in one of the halls of LNGS. The tracking device is made 
of vertical iron walls interleaved with active detectors. The granularity is chosen to 
guarantee an angular resolution of a few degrees for horizontal muons. The overall 
thickness of the iron must be appropriate to effectively separate pions from muons 
at energies greater than 1-2 GeV. These requirements are discussed in details in the 
following sections. 

The signal detection efficiency as well as the rate of pions, muons from tt/K decays 
and muons from semi-muonic charmed hadron decays entering the instrumented surface 
have been computed through a full Monte Carlo simulation based on GEANT 3.21 [To] . 
The target is the LNGS rock, corresponding to a nearly isoscalar target with density 
p = 2.71 g/ cm 3 . Neutrino interactions have been produced by using the event generator 
described in Ref. [TH] and final state particles are propagated through the rock by 
accounting for all physical processes described in GEANT. Particles at the exit of the 
rock are recorded and used for the analysis described in the following. 

3.1 Signal 

An instrumented surface tags v e — > oscillations as an excess of horizontal multi-GeV 
muons plus a small v e — > v T — > fi~X contamination ("silver channel" jTTj) that can be 
safely neglected here. The direction of these muons corresponds approximately to the 
boost direction of the CERN decay ring (~ 3° below the horizon at LNGS) smeared 
by multiple scattering in rock and iron. Their time structure must be consistent with 
the time structure of the circulating beams in the decay ring. A primary muon whose 
vertex is located near the instrumented surface will reach the detector with almost 
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its original energy and in coincidence with other particles belonging to the hadronic 
system, CC interactions occurring deeper in the rock will have a cleaner topology 
due to the screening of the accompanying hadrons but softer muons. Since the energy 
loss of muons in rock is nearly linear with range (~2 MeV g _1 cm 2 ), the target mass 
contributing to the overall event rate grows linearly with the mean muon energy 2 . 
Fig- El (top plot) shows the probability for a muon to exit from the rock with an 
energy greater than 0.5 (full circles), 1 (empty circles) and 2 GeV (empty crosses) as 
a function of the interaction vertex. The muons come from CC interactions in the 
rock and are computed assuming 100% v e — > oscillations at 7 = 2500. Their angle 
distribution with respect to the horizontal direction is shown in Fig. 0] for E^> 2 GeV 
(6 = 90° corresponds to the boost direction of the decay ring). The muon identification 
efficiency as a function of the parent neutrino energy for a rock volume of 40 x 15x15 m 3 
(24.4 kton) is shown in Fig. El (bottom plot). Finally, the number of expected v e — > 
events per year (corresponding to 1.1 x 10 18 18 Ne decays) for 6 13 = 1°,3°,7°, 5 = 0° 
and normal hierarchy as a function of 7 (E^ > 2 GeV) is shown in Fig. 

3.2 Beam related background 

The main sources of background having the same time structure of the signal are 
the punch-through or decayed in flight (DIF) pions from the bulk of v e CC and NC 
interactions in the rock and the semi-muonic decay of charmed particles. Punch- 
through pions are mainly suppressed by the instrumented iron acting as a pion plug 
(iron interaction length Xj = 16.76 cm). The early decays in flight of pions result 
into soft muons which are strongly reduced by the energy cut. Table (0) shows the 
fraction of CC and NC interactions in a 40 x 15 x 15 m 3 rock volume giving a pion 
(muon) that enters the surface with an energy greater than 0.5, 1, 2 GeV. These energy 
cuts correspond to a range in iron of about 2.1, 4.2 and 8.4 interaction lengths. The 
first half of the table refers to 7 = 2500 (u e ) and 7 = 1500 (z/ e ); the second one to 
7 = 4158 (u e ) and 7 = 2488 (p e ). 

The charm production rate has been estimated by using the latest results from 
the CHORUS ^H] experiment and accounting for the different (anti) neutrino energy 
spectrum. For 7 of the order of 2500(z/ e )/1500(z/ e ) the cross sections are 2.3% and 
1.8% for v e and z/ e , respectively. They grow up to 3.5% {v e ) and 2.7% (z/ e ) for 
7 = 4158(^ e )/2488(z/ e ). The semi-muonic branching ratio is about 7%. Note that 
in the highest 7 scenario, the muon spectrum from charm is significantly harder than 
in the region close to the kinematic threshold and the contamination increases sub- 
stantially. This effect is described in Table El where the expected background at 
7 = 2500(f e )/1500(z/ e ) and 7 = 4158(z/ e )/2488(z/ e ) with different momentum cuts are 
shown. It is worth noting that this background would be substantially reduced if the 

2 Similar considerations hold for interactions. However, at a given (anti)neutrino energy, higher 
efficiencies than for CC are expected due to the different y = 1 — E^/E^ dependence of the cross 
section. 
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Figure 3: Probability for a muon (assuming 100% v e — > z/ M oscillations) to exit from 
the rock as a function of the interaction vertex (top panel) and of the neutrino energy 
(bottom panel). The marks represent different energy cuts: larger than 0.5 GeV (full 
circles), 1 GeV (empty circles) and 2 GeV (empty crosses). The /3-beam with 7 = 2500 
is assumed. 
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Figure 4: Muon angular distribution, assuming 100% v e — > oscillations, in z/ M CC 
interactions at the exit of the rock for > 2 GeV. 
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Figure 5: v e — > oscillated events per year versus 7 for #13 = 1°,3°,7 , 5 = 0° and 
normal hierarchy. The rates include the detector efficiency (E^ > 2 GeV) and are 
computed for one 15 x 15 m 2 instrumented surface. 





No cut 


p 7r+( -) > 0.5 GeV 


P7r+(-) > 1 GeV 


p w +(-) > 2 GeV 


u e CC 


1.00(0.62)% 


0.70(0.38)% 


0.41(0.18)% 


0.17(0.07)% 


z/ e NC 


0.76(0.82)% 


0.52(0.54)% 


0.27(0.29)% 


0.11(0.12)% 


u e CC 


0.18(0.41)% 


0.10(0.25)% 


0.05(0.11)% 


0.018(0.024)% 


z/ e NC 


0.33(0.32)% 


0.20(0.20)% 


0.09(0.09)% 


0.02(0.02)% 




No cut 


p 7r +(- ) > 0.5 GeV 


Pk+(-) > 1 GeV 


p n +(-) > 2 GeV 


z/ e CC 


1.59(1.07)% 


1.01(0.71)% 


0.69(0.38)% 


0.35(0.19)% 


z/ e NC 


1.27(1.34)% 


0.90(0.91)% 


0.53(0.53)% 


0.26(0.27)% 


u e CC 


0.36(0.69)% 


0.24(0.46)% 


0.11(0.23)% 


0.04(0.08)% 


z/ e NC 


0.56(0.57)% 


0.38(0.38)% 


0.20(0.19)% 


0.07(0.07)% 



Table 1: Probability for a z/ e , v e CC, NC event generated in a 40 x 15 x 15 m 3 rock 
volume to produce a primary or secondary pion that leaves the rock with a momentum 
larger than a given cut. The first half of the table refers to 7 = 2500 (z/ e ) and 7 = 1500 
(z/ e ); the second half to 7 = 4158 (u e ) and 7 = 2488 (v e ). The first number refers to 
7r + , the second number (in parenthesis) to n~ . 
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No cut 


p„ +( -) > 0.5 GeV 


> 1 GeV 


> 2 GeV 


z/ e CC 


0.030(0.012)% 


0.016(0.004)% 


0.004(0.001)% 


0.001(<0.001)% 


z/ e NC 


0.025(0.018)% 


0.010(0.009)% 


0.002(0.004)% 


0.001(0.001)% 


u e CC 


0.007(0.011)% 


0.003(0.005)% 


0.001(0.001)% 


<0.001((<0.001)% 


z/ e NC 


0.011(0.008)% 


0.004(0.002)% 


0.002(0.001)% 


0.001(<0.001)% 




No cut 


ZV(-) > 0.5 GeV 


Pn+(~) > 1 GeV 


> 2 GeV 


z/ e CC 


0.04(0.02)% 


0.02(0.01)% 


0.008(0.006)% 


0.004(0.003)% 


z/ e NC 


0.03(0.03)% 


0.015(0.014)% 


0.008(0.006)% 


0.002(0.001)% 


u e CC 


0.010(0.014)% 


0.004(0.008)% 


0.001(0.003)% 


0.001(0.001)% 


z/ e NC 


0.018(0.012)% 


0.007(0.004)% 


0.002(0.001)% 


<0.001(0.001)% 



Table 2: Probability for a z/ e , z/ e CC, NC event generated in a 40 x 15 x 15 m 3 rock 
volume to produce a muon that leaves the rock with a momentum larger than a given 
cut. The first half of the table refers to 7 = 2500 (u e ) and 7 = 1500 (z/ e ); the second 
half to 7 = 4158 (z/ e ) and 7 = 2488 (z/ e ). The first number refers to the second 
number (in parenthesis) to 





No cut 


p M > 0.5 GeV 


p M > 1 GeV 


P M > 2 GeV 


Ue CC (7 = 2500) 


6.46 % 


6.33 % 


5.85 % 


2.92 % 


u e CC (7 = 1500) 


5.51 % 


5.51 % 


3.94 % 


0.95 % 




No cut 


> 0.5 GeV 


p M > 1 GeV 


?V > 2 GeV 


z/ e CC (7 = 4158) 


10.63 % 


10.53 % 


10.31 % 


7.63 % 


u e CC (7 = 2488) 


7.15 % 


6.89 % 


5.92 % 


3.44 % 



Table 3: Probability for a semi-muonic decay of charmed v e (z/ e ) CC event generated 
in a 40 x 15 x 15 m 3 rock volume to produce a muon that leaves the rock with a 
momentum larger than a given cut. The first half of the table refers to 7 = 2500 (u e ) 
and 7 = 1500 (P e ); the second half to 7 = 4158 (u e ) and 7 = 2488 (P e )- 



sign of the muon were available, its charge being opposite with respect to the primary 
\x from z/ M CC. For such a soft muon spectrum, a large-surface magnetized iron detector 
with a field of about 1 Tesla would have a rejection factor greater than 99% and, in 
addition, would suppress the punch-through and DIF background by about 50%. 

3.3 Beam unrelated background 

The LNGS underground halls are located at a slant depth that strongly depends on the 
zenithal and azimuthal direction. The minimum depth is ~3200 hg/cm 2 but, due to a 
fortunate conspiracy of the mountain profile, the largest rock depths near the horizon 
are reached in the CERN-to-LNGS direction. The rock coverage in the zenith (9) and 
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Figure 6: Rock coverage expressed in km of water equivalent (km w.e.) in the range 
9 ^ 20° with respect to the beam direction as a function of the zenith and azimuth angle 
(courtesy of the LVD collaboration). The direction pointing to CERN corresponds to 
9 ~ 90° and <p ~ 90° (see text for details). 

azimuth (0) region pointing to CERN is shown in Fig. |HJ Here, 9 = 90° corresponds 
to the horizon; at <p ~ 90° the muons come from the CERN direction and enter the 
instrumented surface from the front. At ~ 270° muons enter the detector from the 
back. The latter can be vetoed if the active detectors provide proper timing. As can be 
inferred from Fig. El the slant depth is greater than 12 km of water-equivalent (km w.e.) 
in most of the region of interest. Here, the muon flux is dominated by atmospheric 
neutrinos and it is of the order of ~ 5 x 10~ 13 cm~ 2 s _1 sr^ 1 [T9~] . A detailed calculation 
has been carried out through a full parametrization of the data from the MACRO 
experiment j20]. Results are shown in Table EJ It is clear that, even at the largest 
angle (9 ~ 40°), the time structure of the beam (10 -4 suppression factor) allows a 
complete rejection of the background. In fact, this result suggests that, opposite to the 
I0W-7 option (CERN to Frejus), the constraint on the bunch length of the beta beam 
(< 14 ns) can be released by at least one order of magnitude. 
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Table 4: Number of underground muons surviving different sets of cuts in 1 year data 
taking with one basic unit target, fi^ stands for the number of cosmic muons entering 
the detector from the front and back side; /if includes just the muons from the front. 
The suppression factor due to the time structure of the beam is not taken into account. 



3.4 Summary of expected rates 

From the above discussion, it is apparent that the sensitivity of the detector under 
consideration will be limited mainly by the beam related background. A severe re- 
quirement on the visible range of the penetrating tracks, corresponding to an energy 
cut of about 2 GeV, will bring the punch-trough contamination at the level of a few 
events per year (suppression factor < 10~ 3 ). The charm contamination is expected to 
limit the sensitivity at the highest 7 (see Table EJ); a further suppression factor of the 
charm background (< 10~ 2 ) from charge reconstruction is available for a magnetized 
detector 3 . In most of the cases, however, the sensitivity is limited by the contamination 
of secondary muons from 7r and K decay in flight. A synopsis of the expected rates 
per year for a 15 x 15 m 2 instrumented surface in the occurrence of the null hypothesis 
(#13 = 0° =$> P[y e — > u^) ~ 4 ) for > 2 GeV and an angle 9 < 40° with respect 
to the nominal beam direction is shown in Table 03 The events per year expected for 
100% u e — > conversion probability are 9.3 x 10 4 (u e at 7 = 2500), 2.0 x 10 4 (P e at 
7 = 1500), 7.9 x 10 5 (u e at 7 = 4158) and 2.1 x 10 5 (p e at 7 = 2488). 

4 Sensitivity 

The sensitivity to the # 13 and 5 parameters is evaluated assuming an instrumented 
surface of 15 x 15 m 2 with a detector having an iron depth greater than 8 interaction 

3 In the few GeV energy range, a magnetized iron detector with B ~ 1 T and, if necessary, precision 
trackers before and after the iron plug can achieve charge misidentification probabilities well below 
10 -2 ^7]. In order to determine the actual efficiency, an optimization of the detector and of the 
pattern recognition algorithm is mandatory. Clearly, this issue is beyond the scope of this paper; in 
the following we assume very conservatively a 99% charge identification efficiency. 
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Table 5: Number of event per year for a 15 x 15 m 2 instrumented surface in the 
occurrence of the null hypothesis (#13 = 0°). B ~ 1 T (B = 0) refers to the detector 
option with (without) magnetic field. 



lengths. Only muons with energy greater than 2 GeV at the entrance of the detector 
are considered. A 5 year data taking with 2.9 x 10 18 6 He and 1.1 x 10 18 18 Ne decays 
per year is assumed 4 . The baseline L corresponds to the CERN to Gran Sasso distance 
(732 km). Results are provided for two high-7 options: 7 = 2500(z/ e )/1500(/> , e ) and 7 = 
4158(z/ e )/2488(z/ e ). The former is dominated by the pion background from u e CC and 
NC. The latter suffers from a significant charm contamination which can be eliminated 
by charge reconstruction (magnetized iron detector option). The sensitivity has been 
computed fixing all parameters but #13 and 5 to their current best values (Am^ = 
7.3 x 10~ 5 eV 2 , sin 2 26 12 = 0.8, \Am 2 31 \ = 2.5 x 10~ 3 eV 2 , sin 2 2 6 23 = 1) and performing 
a two-parameter x 2 fit. The facility is used as a pure counting experiment and the 
information coming from the reconstructed muon spectrum at the entrance of the 
instrumented surface is not exploited. Fig. [7| shows the parameter region excluded at 
90% CL. in the occurrence of the null hypothesis (#13 = 0°) for normal (left plots) 
and inverted hierarchy (right plots) and for 7 = 2500(z/ e )/1500(zv e ) (lower plots) and 
7 = 4158(/y e )/2488(z/ e ) (upper plots). Here, we assume reconstruction of the muon 
charge, use the neutrinos coming from 18 Ne decays (dots), the antineutrinos from 6 He 
(dashed line) and combine the two measurements (solid line). The dot-dashed line 
represents the sensitivity for a detector without magnetic field. Note that, in spite of 
the reduced cross section, the antineutrino flux strongly contributes to the sensitivity 
of the apparatus due to the lower background contamination. As noted above, the 
suppression of matter effect does not allow an unique determination of the neutrino 
hierarchy so that sign [5 ■ sign(Am 2 3 )] remains undetected. The allowed regions at 90% 
CL. for 13 = 1° (red), 3° (black), 7° (green), 5 = 0° (left plot) and 5 = 90° (right plot) 
are shown in Fig. |H] (lower plots) for normal hierarchy, a magnetized detector operating 
at 7 = 2500(z/ e )/1500(z/ e ) and combining both the neutrino and the antineutrino flux. 

4 Note that, in current literature, sometimes results for the low-7 CERN to Frejus option are given 
assuming 10 years of data taking. 
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The corresponding curves for 7 = 4158(z/ e )/2488(z/ e ) are shown in Fig. |H] (upper plots). 
The sensitivity of this apparatus compared with other proposed facilities |3J |^] is 
shown in Fig. El The facilities under consideration are the CERN to Gran Sasso v beam 
(CNGS), JAERI to Kamioka (J-Parc), BNL to Homestake (BNL), the low-7 beta beam 
("Beta beam") and the SPL to Frejus superbeam (SPL). The current CHOOZ limit is 
also shown. The label "This paper" refers to the 7 = 4158(z/ e )/2488(^ e ) magnetized 
detector configuration after 5 years of data taking. 

5 Conclusions 

The v fluxes at large distance produced by j3 decays of boosted radioactive ions has a 
strong quadratic dependence on the Lorentz factor 7 of the ions. For a fixed baseline, 
an increase of 7 that brings the average neutrino energy well above the maximum of 
the oscillation probability does not imply a loss of events since the increase of the 
flux compensates for the smallness of the oscillation probability. Moreover, a net gain 
of events is obtained by the linear rise of the CC cross section and by the increase 
of the detector efficiency The latter effect is particularly valuable for purely passive 
detectors. In this paper we propose to observe v e — > z/ M oscillations as an excess 
of horizontal muons produced in the rock and tracked by a low-mass instrumented 
surface installed in an underground hall of LNGS. This configuration allows a very 
precise determination of #13 and turns out to be competitive with the current low-7 
design based on a gigantic water Cherenkov at Frejus [7j. Similarly to the latter, it has 
limited sensitivity to the sign of Amf 3 . However, opposite to the facilities operating 
at the peak of the oscillation maximum, it shows maximal #13 sensitivity for small 
CP violation (5 ~ 0, ±7r). For 5 = n/2, this complex would significantly constrain the 
9i 3 , 5 parameter space (see Fig. [HJ) but, in general, would not tag explicitly CP violation 
through a lepton/antilepton asymmetry, since the off-peak configuration suppresses 
the CP odd terms. Clearly, the cost of the detector is negligible compared e.g. with 
a Mton water Cherenkov, while most of the acceleration system is either shared with 
other non-neutrino projects (EURISOL) or based on existing CERN machines. 
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Figure 7: 90% CL. excluded region in the occurrence of the null hypothesis 
(#13 — 0°) for normal (left plots) and inverted hierarchy (right plots) and for 7 = 
2500(z/ e )/1500(z> e ) (lower plots) and 7 = 4158(z/ e )/2488(V e ) (upper plots) using the 
neutrinos coming from 18 Ne decays (dots), the antineutrinos from 6 He (dashed line) 
and combining the two measurements (solid line). The dot-dashed line represents the 
sensitivity [y and v combined) for a detector without magnetic field. 
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Figure 8: Allowed regions at 90% C.L. for 13 = 1° (red), 3° (black), 7° (green), 5 = 0° 
(left plots) and 5 = 90° (right plots) for normal hierarchy and a magnetized detector 
operating at 7 = 2500(z/ e )/1500(z/ e ) (lower plots) and 7 = 4158(z/ e )/2488(z/ e ) (upper 
plots). Both the neutrino and the antineutrino fluxes are combined. 
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Figure 9: 90% C.L. excluded region in the occurrence of the null hypothesis (#i 3 = 0°) 
for normal hierarchy and various experimental facilities jHE] (see text for details). 
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